Entanglement dynamics under decoherence: from qubits to qudits 
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We investigate the time evolution of entanglement for bipartite systems of arbitrary dimensions 
under the influence of decoherence. For qubits, we determine the precise entanglement decay rates 
under different system-environment couplings, including finite temperature effects. For qudits, we 
show how to obtain upper bounds for the decay rates and also present exact solutions for various 
classes of states. 
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I. INTRODUCTION 

The production of entangled states and the control 
of their time evolution became a major issue in current 
research in view of the development of quantum infor- 
mation theory, and all possible applications associated 
with it. Besides the formidable experimental advances in 
this direction, there remains a main obstacle which is the 
fragility of entanglement under the unavoidable interac- 
tion with the environment. This coupling of the quantum 
system with its surroundings, and the consequent decay 
of entanglement, motivates important questions such as 
to understand its sources, to identify the characteristic 
timescales, and, possibly, to find ways to circumvent it. 

To devise appropriate strategies for controlling entan- 
gled states under the effect of environment interaction, 
the first step is to acquire a deeper understanding of the 
dynamics of the decoherence processes themselves. De- 
spite the rapidly increasing experimental interest in this 
subject, due to the possibility of monitoring entangle- 
ment dynamics H , most of the theoretical work focused 
on characterizing static properties of entanglement for 
quantum states |^ H, H, |^ - 

Only very recently the question of entanglement decay 
under environment-induced mixing has been addressed, 
for some specific states and environment models, and re- 
stricted to the case of two qubits 0, S ; probably due 
to the lack of a genuine and computable entanglement 
measure for systems larger than that. Nonetheless, new 
techniques for the derivation of bounds 0, 0| of con- 
currence, one possible entanglement measure, recently 
allowed a systematic study of entanglement dynamics 

including multipartite 



and 



for more general states, 
multi-level systems. 

These higher dimensional systems are of great interest 
since they can enlarge the perspectives of efficient ap- 



plications in quantum information and can also be used 
to test fundamental aspects of quantum theory. In fact, 
entangled states of two d-dimensional quantum s yste ms, 
the qudits, can improve measurement resolution jl^ 
and are known to violate local realism more strongly than 
qubits 0, 01 . Moreover they can be used for quan- 
tum computation 0, ItR fl^ [l^ and also in quantum 
cryptography protocols |20ll2lj . which are safer [22, ll^ 
than their qubit counterparts. Despite the importance 
of such entangled qudits, reflected in the intense activity 
on their production and manipulation in different exper- 
imental setups mmillllllillilllOl, their dynamics 
under the influence of environment interaction remained 
unexplored until now. 

The paper is organized as follows. In Section ^] we 
will briefly recall a recently developed approach for cal- 
culating concurrence which allows us to investigate en- 
tanglement between systems of arbitrary dimensions. In 
Section Hill we will present the different models which 
describe the interaction of the system with the environ- 
ment. Section Hvl is devoted to the analysis of the entan- 
glement decay rates under decoherence processes, start- 
ing from the case of bipartite qubits. With the available 
analytical tools, some, still unknown, features of the de- 
coherence dynamics are presented. The section closes 
with the analysis of bipartite qudits. A summary of the 
main results of the paper is presented in the concluding 
Section El 



II. ENTANGLEMENT MEASURE: 
CONCURRENCE 



In order to follow the environment-induced time evo- 
lution of entanglement, one needs a measure which satis- 
factorily deals with mixed states. A commonly used mea- 
sure in the context of two qubits is the concurrence |3l| , 
defined for pure states as 
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where * stands for complex conjugation performed in the 
standard, computational basis. For mixed states it can 
be formulated as 



with 



c{p) = , inf 2^ Pr c(*,) 



Pi > , and p = Vp,|*,)(*, 



(2) 



(3) 



In contrast to most other measures, Eq. |5J can be 
solved algebraically with the well known solution c{p) ~ 
max {Ai — A2 — A3 — A4, 0}, in terms of the square roots, 
Ai , of the decreasingly ordered eigenvalues of the matrix 

P(cr.y (Kl ay)p*{ay (g) (Ty). 

For higher dimensional systems, we will use the genen- 
eralization of concurrence given in 32] that coincides - 
though not obviously - with the original concurrence 
if restricted to two-level systems. For practical purposes, 
it is convenient to express concurrence [ssL Is^ in terms 
of an operator A acting on the product space Ti^Ti. oi 
two copies of the system, as 



The operator 



(4) 



(5) 



is the projector onto the space spanned by the states \xa) 
that are anti-symmetric with respect to the exchange of 
the copies of either TLi or TL2, i-e. 



\Xa) = i\ikil) - \ilik)) «> {\jrnjn) - \jn jm)) 



(6) 



The states and {|jm)} form, respectively, arbitrary 

local bases of 7ii and H2, and a is a label for the multi- 
index [fc, Z, TO, n\. 

To extend this construction for the case of mixed 
states, one should substitute Eq. Q in the convex-roof, 
Eq. ^ , which can be written, in terms of subnormalized 
states \ipi) = y^il^i), as 

C(p) = inf V J{yJ^\®{^^\ A \4'^)®\^P^). (7) 

{i>i} ~' * 

Unfortunately, apart from the two-level case where the 
operator A can be written in terms of a single vector 
1;,^) = (|01) - |10)) (E) (|01) - |10)), no exact solutions for 
Eqs. (O and ^ are known and, hence, one has to rely on 
numerical efforts to calculate the concurrence. Also note 
that numerical solutions of this optimization procedure 
define only an upper bound for concurrence, since there 
is no a priori information available on whether the global 
minimum or just a local one has been reached. 

To circumvent this problem, we will use suitable ap- 
proximations which provide lower bounds of concur- 
rence 0, ll(ll | and not only allow for an efficient numer- 
ical approach, but also, in some cases, for exact alge- 
braic solutions. Starting from a decomposition of the 



density matrix in terms of pure (subnormalized) states 
p — J2i \4'i){4'i\ ctud from the vectors \xa), one can de- 
fine a set of matrices T", given by 



T^k^ix^f^l^k) 



(8) 



These are connected to the previously defined operator 
A through the tensor 



It was shown in j^] that concurrence is bounded from 
below by 



c(/o) > max < Si 



i>l ) 



(10) 



where the Si are the decreasingly ordered singular values 
of 



with^^ 



1 



(11) 



that still depend on the choice of the complex parameters 
Za- Although any choice provides a lower bound, one 
might still wish to carry out an optimization over Z^, 
though now on a much smaller parameter space and with 
simpler constraints than in Eq. I?)). Moreover, also each 
matrix T" already provides a lower bound, which can be 
calculated algebraically. 

Finally, let us describe an experimentally motivated 
approach to calculate lower bounds of concurrence, the 
quasi-pure approximation p|. Although environmental 
influences cannot be avoided completely, under typical 
experimental conditions one deals with states which are, 
at least initially, quasi-pure. This is, they have a sin- 
gle eigenvalue fj,i which is much larger than all the oth- 
ers, and an approximation based on this condition can 
be developed. Indeed, using the spectral decomposition 
p — J2i of the density matrix and the previ- 

ously defined subnormalized states, one can see that the 
elements of A defined in Eq. © are proportional to the 
square roots of the eigenvalues p,i : 



(12) 



This proportionality and the assumption that pi ^ pi 
allows to order the elements of A in powers of the square 
roots of Pi. Keeping only the leading order terms, i.e., 
A\{ at zero order, ^j^, ^Jj, A}j\ and A^j at first order, 
and so on, we can approximate A by 



A Im 



jk 



with 



r(qp) 
^jk 



■^jk 



(13) 
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Thus the quasi-pure concurrence can be written as 



cip) 



{p) = max < 5i — Si 



(15) 



i>l 



where the Si are the singular values of the matrix T^'ip) 
defined in Eq. ((Tl|l . 

With these tools at hand, entanglement in higher di- 
mensional systems can be explored in a computationally 
manageable way. We will use them to monitor the entan- 
glement dynamics under different sources of decoherence. 



III. ENVIRONMENT MODELS 

We start out with the assumption that each subsys- 
tem interacts only, and independently, with its local en- 
vironment and, therefore, any initially entangled state 
will evolve, in some cases asymptotically, into a separa- 
ble one. The dynamics of this decoherence process can 
be described by the master equation 



(16) 



where p is the reduced density operator of the system. 
The interactions of each component of the system with 
the environment are assumed to be of the same form, rep- 
resented by the Lindblad operator C. Under the assump- 
tion of complete positivity and Markovian dynamics, the 
action of £ on p reads 



(17) 



The operators Li and the rates Ti define, respectively, the 
specific type of the system-environment coupling, and its 
strength. Various physical situations may arise when a 
system is coupled to a reservoir: dissipation can take 
place, noise can be added to the system, or simply loss 
of phase coherence can occur. All these processes corre- 
spond to a suitable choice of the operators Li in Eq. H17|) 
and lead to entanglement decay, with a rate which is at 
the focus of this paper. 

For two-level systems, the operators Li can be written 
in terms of the Pauli matrices. In the case of an interac- 
tion with a thermal bath the Lindblad operator reads 

r(n + 1) 

Cp = ^ (2 a^pa+ — a+a^p — pa+a^) + 

— (2(T+p(T_ - cr_(T+p- pcr_(T+) . (18) 

In this equation, the first and the second term on the 
right hand side describe, respectively, decay and excita- 
tion processes (mediated by the two-level excitation and 
deexcitation operators fT+ and cr_), with rates which de- 
pend on the temperature, here parametrized by n, the 
average thermal excitation of the reservoir. In the zero 



temperature limit, n = 0, only the spontaneous decay 
term survives, leading to a purely dissipative process, 
which corresponds to the fundamental limiting factor for 
the coherent evolution of atomic qubits. Noisy dynam- 
ics corresponds to the infinite temperature limit, where 
n — !■ oo, and, simultaneously, F — s- 0, so that Tn = F 
is constant. In this case, decay and excitation occur at 
exactly the same rate, and the noise induced by the tran- 
sitions between the two levels brings the system to a sta- 
tionary, maximally mixed state. Finally, a purely dephas- 
ing reservoir is obtained by choosing Li — d — a+U-, 
leading to the master equation 



dp 
'dt 



(2 a+(J^pa+a^ — a+a^p — p(y+cr- 



(19) 



In this case, only the off-diagonal elements of the density 
matrix decay, and phase coherence is lost. This dephas- 
ing mechanism is known to be an important source of 
decoherence in ion traps |l35t. .36j and in quantum dots 
experiments |37ll3ll39|. 

If we consider a system of qudits, the operators Li in 
Eq. IjlTI) cannot be written anymore in terms of Pauli 
matrices. To describe the effect of the environment in 
such systems, we will consider the qudits as two bosonic 
modes with truncated bases of length d. The terms n and 
TO in the general (pure) state \^) = Ylt^m=o'4'nm\nm) 
then indicate the occupation number in each of these 
modes. Decay and excitation processes now correspond 
to the action of annihilation and creation operators a and 
a\ analogous to the (t_ and (t_|_ for qubits. With this 
convention, the thermal and dephasing environments are 
described, respectively, by 



Cp 



Tin + 1) + + 4. N 
^ (2 apa' — a' ap — pa' aj + 

Tn 

{2 a} pa ~ aa) p — paa}^ , (20) 



and 



Cp — — {2a'^ apa^ a ~- a^ ap ~ pa^ a^ . (21) 

Zero and infinite temperature limits are obtained as in 
the qubit case, and represent, to a very good approxi- 
mation, the decoherence processes of photons in high-Q 
caviti es l iol . and of motional degrees of freedom in ion 
traps |4l|. 



IV. ENTANGLEMENT DYNAMICS UNDER 
DECOHERENCE 

A. Two qubits case 

We shall begin our analysis with the case of two qubits, 
since it allows for analytical solutions, and can provide 
some intuition of the processes at work in higher dimen- 
sional systems. To provide a more complete and system- 
atic description of the timescales of entanglement decay 
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under environment interaction, we will recall some known 
results for the case of dephasing and zero temperature 
reservoirs, and furthermore derive exact solutions when 
finite temperature effects are taken into account. 

Let us start with initially prepared Bell states l^**) — 
(|01) ± |10)) /V2, and = (|00) ± |11)) /V2. Their 

concurrencies as a function of time are given by 

c(*±,i) =c($±,i) = e-", (22) 

for the dephasing environment, and by 

c(*±,i) = e-", c($±,t) = e-2r*, (23) 

for the zero temperature case. The first important ob- 
servation is the accelerated (by a factor of two) decay of 
concurrence for the as compared to the states, 
under the influence of a zero temperature environment. 
This can be understood from the timescales involved in 
the corresponding solution for the density matrix: while 
for each term |01) and |10) corresponds to a single 
particle decay, leading to a timescale e~'"*, we have the 
term |11) in such that both particles can undergo an 
environment induced transition to the ground state, thus 
introducing a faster, e"^'"*, decay. 

Also for the finite temperature case can an explicit 
solution be calculated for initial Bell states: it reads 



c{t) — max {cT(i), 0} 



(24) 



where the function exit), for ^E"* and states, is given, 
respectively, by 



CT(*±,t) = /?- 

and 



2(1 - 13) ^{n^ + n)2 (/3 + if + ^{ffi + n) 



{2n + lf 



(25) 



[2n{n +l) + l](3^ - (3- 2n{n + 1) 
{2n + iy 



(26) 

with [3 ^ e r(2ri+i)t q;^j^ggg expressions, which reveal the 
precise role of temperature on the decay of entanglement 
of Bell states, show that, distinct from the cases of de- 
phasing and zero temperature environments, concurrence 
does not decay according to a simple mono-exponential 
law, but rather exhibits various timescales. To apprehend 
the essential features of temperature effects on entangle- 
ment decay, it is useful to focus on the short time limit 
of these equations Expanding Eqs. (|25|l and H26() to 
first order in i, one finds that the short time concurrence 
decay is given by 

c(^'±,i)~l - (2n+l + 2v/n(rl + l))n, (27) 



and 



This shows that, at the beginning of the evolution, a fi- 
nite temperature reservoir increases the concurrence de- 
cay rate, as compared to the zero temperature case, by 
a factor of {2n + 1 + 2^n{n + 1)^ and of (2n + 1), for 
^1'='= and states, respectively. 

From Eqs. H25|l and H26fl one can also derive the infinite 
temperature limit (n — > oo, F ^ 0, with Tn = F = 
const.) of the thermal bath. In this case, the concurrence 
becomes 



c(vE'±,i) = c(<i>±,i) 



-4rt 



(29) 

Note that, again, entanglement dynamics involves differ- 
ent timescales. However, a single exponential captures 
the basic behavior of entanglement decay, since Eq. H29() 
can be well fitted by a function in the form ae~^* + S, 
with negative offset S. Consequently, separability, i.e. 
c{t) = 0, is reached at finite times for the infinite tem- 
perature environment, as well as for any n > 0. This is 
in contrast to the above dephasing and zero temperature 
environments, where separability is reached only asymp- 
totically, and follows from the long-time limit for exit) 
in Eq. ^ 



lim CT{t) 

t — >oo 



2n{n + 1) 
' {2n+lf 



(30) 



c($=^,t) ~ 1 - 2{2n + l)Tt. 



(28) 



which is non-positive for all n > 0. The exact expres- 
sions for this separability time can be easily obtained 
from the condition c{t) = 0, and are of the form tg^p = 
In [/(n)]/F(2n -|- 1), where the function /(n) depends on 
the initial state Vt* or Observe, however, that this is 
not true for arbitrary initial states. Indeed, some initially 
mixed states Q as well as some pure non-maximally en- 
tangled states, e.g., = i|00) + + reach 
separability on finite timescales even for a zero tempera- 
ture environment. 



B. Two qudits case 

As the dimensions of the subsystems increase, not only 
a general analytical solution is unknown, but also the nu- 
merical approach to obtain reliable estimates of concur- 
rence rapidly turns into a very demanding task. How- 
ever, it is exactly in this situation where the strength 
of the tools derived in Section ^ for calculating lower 
bounds of concurrence, becomes manifest. We will de- 
vote the remainder of this paper to a systematic analysis 
of entanglement dynamics in bipartite qudits, using these 
tools. In particular, we will extract bounds for the decay 
rates of entanglement, and also show how one can infer 
analytical solutions even for arbitrary dimensions, from 
the knowledge of the dynamics of the system. 
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1. Dephasing environment: an exactly solvable example 

We begin our analysis with the dephasing dynamics 
described by Eq. H21|l . which, as discussed in Section ITTll 
induces a decay of the off-diagonal elements of the density 
matrix without changing the diagonal ones. Besides its 
importance in the description of different physical situa- 
tions, as already mentioned before, the dephasing model 
is very instructive since it allows for obtaining analytical 
solutions for the concurrence of some classes of states. 
Let us first assume that the initial pure state is of the 
form 1-0) = ajmi TO2) + b\nin2), what comprises two- 
level Bell states as special cases. In this situation, one 
can write the solution of Eq. H21|l as 

p{t) = \a\^\mxm2){mim2\ + |6|^|?^.l?^2)(f^lf^2| 
+e~'^*'ah*\mim,2){nin2\ + e^^*a*h\nin2){mim2\, (31) 

with 

l^^[{mi-nif + [m2-n2f]. (32) 

Note that there is just one decaying element (and its 
conjugate) in the density matrix evolution and, therefore, 
one expects that the entanglement decay rates should 
follow the same behavior. This is indeed the case, as we 
will show using the techniques described in Section ^1 
First, one can easily check that the mixed state (|31|l can 
be decomposed into contributions of two pure states, 

= y/p{a\mim2) + b\nin2)) , (33) 

and 

102) = V"! -p{a\mim2) - b\nin2)) , (34) 

withp= (l + e-T*)/2, and |0i>(0i| + |02)(02|- From 
the structure of this decomposition, one can see that, 
amongaU Ixq), only Ix) = (|mirii) - |nimi))(g) (|to2"-2) - 
|n2TO2)) gives a non-zero contribution for T" in Eq. ||SJ|. 
In this case, the concurrence can be expressed in terms 
of a single matrix T, which reads 

{ 2pab \ 
2{p~l)ab J ■ 

Thus the formal structure of two-level systems is re- 
trieved, and the lower bound, Eqs. ((Tn|l and ((TTl) . is exact. 
From the singular values of the matrix T, one readily ob- 
tains 

C{t) = 2\ab\ e(-rt/2)[(mi-ni)^ + (m2-n2)^] ^ (35) 

which is equal to the sum of the absolute values of the 
two non- vanishing off-diagonal elements of p(t). 

2. Zero temperature environment: exact solutions, and 

estimates 

Although there are no a priori criteria from which to 
infer the existence of exact solutions like Eq. H35|) in gen- 
eral, this is still possible for a large variety of physical 



situations. For a zero temperature environment, for ex- 
ample, exact solutions can be derived based on a given 
form of the decomposition of p. Moreover, even if such 
decomposition does not exist, it is still possible to obtain 
estimates on the decay rates of concurrence with the tools 
described in Section im 

To illustrate this, we start with states of the form jV^) = 
a\Om) 4- b\mO) and consider, for the derivation of the 
exact solution, the case of three-level systems and a — 
b — l/\/2. The decomposition of p{t) in terms of its 
eigenstates (with non-vanishing eigenvalues) reads 

|</.i)=e-i^*(|02)-f |20))/V2, (36) 

102) = ^e-r*(l-e-r*)|01) , (37) 

103) = y^e-r'(l-e-r*)|10) , (38) 

and 

104) = (1 - e-r*) |00) . (39) 

All these states, with the exception of |0i), are separable. 
Therefore, from Eq. (|2Jl, one can see that the concurrence 
of |0i) is an upper bound of c{p). Moreover, this term 
is also the only one that gives a non-vanishing contri- 
bution to the matrix r('ip), Eq. 114|) . in the quasi-pure 
approximation. Hence, lower and upper bounds coin- 
cide and the solution c(t) = e~^'"* is exact. For the state 
10) — a\Om) +b\mO) the demonstration is analogous, and 
the concurrence is given by 

c(t) = 2|a6| e-"^*. (40) 

From the time-dependent matrix elements of p, which can 
be obtained from Eq. (|20|l with n = 0, one can see that, 
as in Eq. H35() . concurrence is simply given by the sum 
of the absolute values of the only two non-vanishing off- 
diagonal elements of p. However, we should emphasize 
that this only happens in special cases alike the present 
one and, in general, there is no simple relation between 
the decay of off-diagonal elements of p and concurrence. 

An example where such a simple relation does not ex- 
ist, even though an exact solution is possible, is the zero 
temperature dynamics of initial states of the form |0) = 
a|00) -|- b\mm). The time-dependent density matrix can 
be decomposed as p = ^+r], with rj diagonal in the {|« j)} 
basis (and thus separable), and ^ — J2{^j\p\Pl)\^j){Pl\ 
with i, j, p, q £ {0,to}. From this argument and from 
Eq.(|2Il it follows that the concurrence of ^ is an upper 
bound of c{p). Furthermore, ^ is a two-qubit-like matrix 
and its concurrence can be readily obtained as 

c{t) = 2e-"i^* {\ab\ - (1 - e-i^*)" \b\^) . (41) 

This result coincides with the quasi-pure approximation 
of c(^), calculated using \x) = (|0r7i) — |toO)) ® (|0m) — 
|mO)) in Eqs.ljHJ, lEJl, ((T^ . and hence is exact. 
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It is clear from Eq. (|41() that the entanglement decay, 
in this case, encompasses different timescales for d > 2, 
while the only non-zero off-diagonal elements, {00\p\mm) 
and (mm|p|00), decay as e~™'"*. Note also that, in 
the limit of large m, the concurrence tends to c(t) — 
2\ab\e~™^* and then coincides with the result of Eq. (|^ . 
Thus, in this limit, the state = a\00) + b\mm) decays 
as the state = a\Om) + b\mO) . This shows that the ro- 
bustness of the Bell states as compared to the in 
terms of decay rates in the two qubit case (see Eq. 
fades away for larger dimensions. 

In the two previous examples, exact solutions were pos- 
sible because p could be decomposed into a separable 
part, and a non-separable two- level- like one whose upper 
and lower bounds coincide. Although, at first sight, one 
might judge these conditions to be rather restrictive, the 
above considerations can be also applied to other envi- 
ronment models and qubit-like initial states. However, 
what happens when such decomposition does not exist? 
How much information on the decay rates can we extract 
with the available tools? 

To help to answer these questions let us focus on 
the situation shown in Fig. ^ where an analytical 
solution could not be found. For the initial state 
\iP) = (|12) -I- |21)) /v/2, upper © and optimized lower 
bounds IjlUI as well as quasi-pure calculations (|15|l . 
are presented in the upper panel of the figure, as cir- 
cles, crosses, and squares, respectively. The solid lines 
represent the quasi-pure approximations for initial states 
of the form = {\lm) + |ml)) /\/2, with m = d — 1 
and, from top to bottom, d = 4 to d = 7. Note that 
the quasi-pure approximations are calculated only until 
times when the initially largest eigenvalue Ai coincides 
with the second one (see the lower panel of Fig. ^ for the 
case d = 3). The quasi-pure solutions solutions have a 
simple form e^(™+i)r*_ Por d = 3, we can go further and 
use also the numerically calculated upper bound, which 
in this case is found (through curve fitting) to behave as 
e~^^*, to confine the actual value of the decay rate to 
the interval 2r < 7 < SF. Furthermore, we can state 
that for arbitrary dimensions the decay rate cannot ex- 
ceed the value given by quasi-pure approximation, which, 
in the present case, leads to 7 < dT. This is a notice- 
able asset of our analysis. Even if the high dimension of 
the system, together with the complex structure of the 
states, do not allow to determine the precise decay rates 
of entanglement in the general case, we can extract useful 
information on their bounds. 



3. Finite temperature environment 

Finally, let us briefly present a situation where the abil- 
ity to handle higher dimensional systems is essential: the 
finite temperature reservoir. Any initial state will evolve 
asymptotically into a thermal state, which can be de- 
scribed by a finite number of levels - which, however, 
increases with temperature. 




0.25 0.5 0.75 1 

Ft 



FIG. 1; Top panel: optimized upper (circles) and lower 
(crosses) bounds of concurrence, as well as quasi-pure approx- 
imation (squares), for the initial state \il>) = (|12) -I- j21)) /\/2 
coupled to a zero temperature reservoir. The solid lines 
show the quasi-pure approximations for initial states \ip} = 
(jlm) -I- jml)) /\/2, with m — and, from top to bottom, 

ci = 4 to d = 7. Quasi-pure solutions behave as e"'-™"'"'^-"" ' 
and the upper bound for d = 3 is correctly fitted by e"^""*. 
Bottom panel: eigenvalues of the density matrix p{t) as a 
function of time for d — 3. Quasi-pure results coincide with 
the optimized lower bound and are only calculated while Ai 
remains the largest eigenvalue. 



The decoherence effect of this environment on qu- 
dits is illustrated in Figure 13 for an initial state \tfj) = 
(|01) + \10))/V2. The solid, dashed and long-dashed lines 
represent, respectively, the two qubit solution, Eq. if^ . 
for fi = 0.1, n = 0.2, and the infinite temperature limit. 
The corresponding quasi-pure solutions for the same ini- 
tial state evolving under Eq. (|20|l are given by circles, 
squares and crosses. The bold solid curve indicates, for 
comparison, the zero temperature solution. The dimen- 
sion used in the simulations (d = 8) was chosen to be 
large enough to describe the dynamics of the system for 
n = 0.1 and n = 0.2, for the times shown in the figure. 
Note that for the infinite temperature case actually an 
infinite basis is needed and that our truncated basis set 
allows for a correct description of the dynamics only for 
a short period of time (F t « 0.06), as long as the dynam- 
ics do not populate levels at the boundary of the Hilbert 
space. 

The most noticeable effect here is that the expected en- 
hancement of the decay rates with temperature is more 
pronounced for qudits than for qubits. As a matter of 
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FIG. 2: Finite temperature effects on the decay of entangle- 
ment for an initial Bell state Solid, dashed and long 
dashed lines are, respectively, the two-level case solutions of 
Eq. CHJi given by Eq. for n = 0.1, n = 0.2, and the 
infinite temperature limit. Circles, squares and crosses are 
the corresponding quasi-pure results for the qudits' dynamics 
given by Eq. ^ using a truncated basis with d = 8. For 
comparison, also the solution for a zero temperature environ- 
ment is shown (bold line). The expected enhancement of the 
decay rates with increasing temperature is more pronounced 
in the higher dimensional case. 



fact, this is not too surprising since the dynamics of a 
finite temperature bath, given by Eq. (|20|l . induces tran- 
sitions to many different levels, thus speeding up entan- 
glement decay. 



V. CONCLUSIONS 

In this paper we have contributed to the understand- 
ing of the dynamics of entanglement under environment 
coupling for systems of two qudits. Although no analyti- 
cal solutions are available in general, we have shown that 
they can be derived for different classes of initial states 
and dynamics. Moreover, we have shown how suitable 
decompositions of the density matrix, together with tech- 
niques for calculating lower bounds of concurrence, can 
be useful to decide on the existence of exact solutions 
and to derive them for concurrence or its lower bounds. 
For more general cases in high dimensional systems, we 
were able to derive upper bounds for the decay rates of 
entanglement. 

Future work will have to focus on error estimates for 
the employed approximations for concurrence, in higher 
dimensions. 
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